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The phase transition between the nuclear matter and the quark matter is examined. The relativis- 
tic mean field theory(RMF) is consider with interacting nucleons and mesons using TM1 parameter 
set for the nuclear matter equations of state. It is found that the trasition point depends on cou- 
pling constant a s and bag pressure. From the study of the structure of a hybrid neutron star, it 
is observed that the star contains quark matter in the interior and neutron matter on the outer 
perifery. 
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\ It is widely believed that nuclear matter undergoes a phase transition to quark matter at high densities and/or high 
temperatures. The high temperature limit is expected to have interesting consequences in heavy ion collision and/or 
in cosmology, whereas high baryon density behaviour is important for the study of neutron stars. 

It is expected that Quantum Chromodynamics (QCD) as the fundamental theory of the strong interaction should 
explain possible modifications of hadron properties in the nuclear medium. However, typical nuclear phenomena 
at intermediate and low energies cannot be analytically derived from QCD although one hopes that QCD will be 
solved numerically on the lattice in near future. Meanwhile we are left with the construction of phenomcnological 
models in order to try to describe nuclear phenomena and its bulk properties. Walecka and others |l],|2j used a kind of 
relativistic scalar-vector theory to describe the nucleon-nucleon properties of nuclear matter as well as of properties 
^—1 ' of finite nuclei. Some of the drawbacks of the original Walecka model are that the effective nucleon mass obtained at 
high densities is too small and its incompressibility at the energy density saturation is too large. To eliminate these 
0^ 1 difficulties Boguta and Bodmer ||] modified the original model by introducing selfcoupling terms to the scalar field. 

The inclusion of nonlinear terms to the scalar field surprisingly improved the results of nuclear matter as well as of 
finite nuclei. However, in most of the successful parameter sets, the last term of the self coupling constant is found to 
be negative Q . This negative value of the last term gives an unphysical situation at high density, which is essential 
for a further modification of the model. This modification is done by Bodmer || by introducing a quartic nonlinear 
term to the vector potential to study the equation of state. Later on this suggestion was considered to study the 
properties of finite nuclei [|| , which gives the nuclear matter and finite nuclei properties excellently well (?j . 

Many of the existing parametrisation is unable to reproduce the properties of finite nuclei, nuclear matter and that 
of the physical properties of accreting steller objects, like neutron star etc. Recently a detail calculation has been 
done with different models to study the properties of neutron stars || . Glendenning || has studied the properties 
of neutron star in the framework of nuclear relativistic field theory. This parameter set is not applicable directly to 
finite nuclei. Similarly many parameter sets which explain the properties of accreting matter is unable to explain the 
properties of finite nuclei as well as of normal nuclear matter properties. On the other hand, the parameter set which 
is able to explain the properties of finite nuclei and normal nuclear matter properties, failed to explain the properties 
of accreting matter. In this work, our aim is to see the applicability of the improved parametrisation of Sugahara and 
Toki. Q to nuclear matter which explains well the properties of finite, including superheavy nuclei, infinite nuclear 
matter and the properties of negative energy bound states at normal as well as at high densities Jl0| . 

In a similar study [jil"| , the equation of state for neutron matter is obtained in a nonperturbative method with pion 
dressing of neutron matter, an analysis similar to that of symmetric nuclear matter. The quark matter sector was 
treated perturbatively with bag constant Bi = 148 GeV. Stable solution for such a quark-neutron hybrid star was 
obtained with Chandrasekhar limit as 1.58A/q and radius around 10 km. However in the present calculation we have 
considered the effect of a — lu mesons with the nonlinear interactions and have observed the increase in size and mass 
of such stars with mass about 2.3Mq and radius 13.5 km. 

In a further study, Sugahara and Toki (lj] have taken A — u> tensor coupling and found a heavier critical mass 
of neutron star beyond observational boarder but without tensor coupling they have shown that their result does 
not agree with observational result indicating that the tensor coupling is indespensable for meeting observational 
requirement. 
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We consider here relativistic mean field theory with interacting nucleons and mesons, using a nonlinear version 
in both a and lo mesons for the nuclear matter equation of states. Quark matter is treated perturbatively for high 
densities at short distances |l3| . A first order phase transition between nuclear matter and quark matter seems to be 
indicated. Solutions of the Tolman - Oppcnhcimcr - Volkoff (TOV) equations yield a hybrid star having a quark core 
with a crust of neutron matter. 

The paper is organized as follows. In Sec. II, we present a brief theory for nuclear matter (neutron matter) equation 
of state. The quark matter equation of state is discussed in Sec. III. In section IV, we discuss the structure of hybrid 
neutron star. A summary and concluding remarks are given in Sec. V. 



II. NUCLEAR MATTER EQUATION OF STATE 

We start with the effective Lagrangian density for a nucleon-meson many-body system for nuclear matter. In this 
Lagrangian we have considerd only the interction of nucleons with a , lo and p mesons. 
The Lagrangian is given as jj]]4||G|] 

C = ^(iffy - M)ip + ^d»ad»o - ^mla 2 + ^g 2 a 3 + ^g 3 a 4 - g s ^a 

- ~W U % V + \mlcj^ + \c a (u^f ~ g^^Lo, - ^R a ^R a ^ 

The fields for the a,u) and p-mesons are denoted by a, lo^ and Rp respectively and ip is the Dirac spinor for the nucleon. 
Here g Sl g w , g p are the coupling constants for a, lo and p- mesons and 172, 93 and C3 are self coupling constants. M 
is the mass of the nucleon and m a , and m p masses of the a, lo and p-mesons respectively. The contribution of 
p-mesons to neutron matter is essential and has effect on the formation of hybrid stars . 

In Ref. it has been shown that RMF approach is successful to describe the result of Relativistic Dirac Bruckner 
Hartree-Fock (RDBHF) calculations in nuclear matter. It is found that although the RMF model with scalar self- 
interactions is able to describe effectively the binding energy of nuclear matter as well as the bulk properties of finite 
nuclei, this is not followed by a proper description of the effective nucleon mass and a time-like component of the 
vector self-energy. This is caused mainly by a too restrictive treatment of lo— meson term in the RMF approach, which 
does not take into account the density dependence produced by the relativistic Dirac-Brueckner approximations to all 
mesons involved in the theory. The study of finite nuclei and nuclear matter of Sugahara and Toki H shows that the 
vector potential of the RMF theory increases linearly with density and gets stronger, while RDBHF bends down with 
density. The scalar potential of the RMF theory seems to be overestimating the RDBHF results at high density in 
order to compensate for the strong repulsion in the vector channel. This is the reason for providing the wrong sign in 
ex 4 self-coupling constant in most of the successful parameter sets. Thus, Sugahara and Toki introduced a nonlinear 
term (lo^lo^) 2 into the lo vector meson potential to study the properties of finite nuclei. 

In the mean field approximation the meson field operators are replaced by their expectation values. We also consider 
the isotropic system at rest. The equation of motions for meson and nucleon fields are 

m 2 a a = -g s p s - g-io 2 - g 3 <r 3 (2.2a) 
ttiIloo = g^PB - c 3 lOq (2.2b) 

™ 2 P Rl = \g P < ^7°^ > (2.2c) 
m 2 p Ro3 = g P P03 (2. 2d) 

/>l)3 = ± < ^f T ^ > (2.2c) 



where 

1 



and 

{-%&-y + pM*)il> = {E-g u u>o)1) (2.2f) 

In the above we use the effective nucleon mass M* = M + g a a. The source terms for scalar and vector fields are the 
scalar density p s =< ijrip > and the vector(baryon) density ps =< > respectively. Using the standard positive 
energy solutions of the Dirac equation, we obtain 
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Here we assume that nuclear matter consists of filling nucleon levels up to the Fermi momentum kf and 7 is the 
spin-isospin degeneracy factor (7 = 4 for nuclear matter and 7 = 2 for neutro n mat ter). The effective nucle on m ass 
M* has to be determined self-consistently at each density by solving equation ( 2.2a ) for the scalar field and (|2.3| ) for 
the scalar density. 

The energy density of the nuclear matter in the mean field approach is given by 



e = e N + e c 

Here is the energy density of nucleons of mass M* 
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The energy density e a is the sigma meson interaction term which may be written as 

-^22^ 3 4 

The third term e u is the omega-meson interaction term which is given by 

1 2 .2 1 1 
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Similarly the pressure for the nuclear matter is given by 

P = P N + P a + P u + P p 



(2.6d) 



(2.7) 



Pn = - 



1 7 



3 2tt 2 

iJL 

3 2tt 2 



Vfc 2 + M* 



4 / 



A: 2 



Q J 



k 2 + M* 2 + -M* 4 ln 



kf + Jkj + M* 2 



M* 



(2- 



1 9 9 1 3 1 4 

Pa = —^rn a a ~ 3-9 2Cr ~ 4^3^ 
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Here we use TM1 parameter set. The values of the parameter set are || M = 938.0, m a = 511.198, m w = 783.0,m p 
= 770, g s = 10.0289, g w = 12.6139,5 P = 4.6322 g 2 = -7.2325 fm" 1 ,^ = 0.6183 and c 3 = 71.3075 where masses 
are in MeV. The corresponding nuclear matter properties obtained from the parameter set are po — 0.145 fm -3 , 
E/A = -16.3 MeV, K = 281 MeV, M*/M = 0.634 and a asy = 36.9 MeV, where p , a asy , K, and E/A are the density, 
asymmetric parameter, compressibility modulus and the binding energy per particle, respectively. Sugahara and Toki 
H show that the TM1 parameter set gives more closer results with the RDBHF than the NL1 and NL-SH parameter 
sets. Also, we know earlier g] that the linear set of Horowitz and Serot gives stiff equation of states and predicts 
a too high value of compressibility modulus of about 560 MeV whereas the empirical value is 210 ± 30 MeV [ pr[ . 
The value obtained by TM1 parameter set is more convincing. In figure 1, we have shown the behaviour of pressure 
against density in TM1 and Nl-SH parameter set. 
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FIG. 1. Pressure (P) versus density (p) of nuclear matter is shown with TM1 parameter set. 



III. QUARK MATTER EQUATION OF STATE AND PHASE TRANSITION 



Existence of quark matter in the core of neutron stars/pulsars is an exciting possibility Jl5| . Densities of these 
stars are expected to be high enough to force the hadron constituents or nucleons to overlap thereby yielding quark 
matter. Since the distance involved is small, perturbative QCD is used to derive quark matter equation of state. We 
take the quark matter equation of state as in Refs. jl6 17 1 in which u,d and s quark degrees of freedom are included 
in addition to electrons. Here we set the electron, up and down quark masses to zero Q and the strange quark mass 
is taken to be 180 MeV. In chemical equilibrium fj,^ = fj, s = fj, u + fj, e . In terms of baryon and electric charge chemical 
potentials \jlb and he, one has 
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(3.1) 



Confinement is simulated by a bag constant B. The electron pressure is [HQ 



The pressure contributed by the quarks is computed to order a s = g 2 /47r where g is the QCD coupling constant 



(3.2) 
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The pressure for quark flavor f, with f=u,d or s is jl^]lrj , |l7jl 
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The Fermi momentum is kf = — m 2 ) 1 / 2 . The total pressure, including the bag constant B is 

P = P e + J2Pf-B. (3.4) 
/ 

There are only two independent chemical potentials [Ib and /ie- /*b is adjusted so that the matter is electrically 
neutral, i.e. dP/dfiE = 0. The baryon number density is given by p — dP/dfiB- 

We now consider the scenario of phase transition from nuclear matter to quark matter. As usual, the phase 
boundary of the coexistence region between the nuclear and quark phase is determined by the Gibbs criteria. The 
critical pressure and critical chemical potential are determined by the condition 

Pnm(HB) = Pqrn(HB)- (3.5) 

We take a s = 0.5, 0.6 and the bag constant B = (150 MeV) 4 , (155 MeV) 4 , which is a reasonable value to calculate 
pressure in the quark sector. Schaab et al. ]l8| used this value to be i? 1 / 4 = 145 MeV. However in a calculations 
of Glendinning the bag pressure was taken as i? 1 / 4 = 180 MeV. In this calculations |l9| the transition was 
determined for the above bag constant which places the energy per baryon of strange quark matter 1100 MeV, well 
above the energy per nucleon in infinite nuclear matter as well as the most stable nucleus, 56 Fe (E/Am 930 MeV). In 
figure 2,we have plotted pressure versus chemical potential for nuclear matter and quark matter. The solid line is 
shown for nuclear matter,with TM1 parameter set. The dash and the dotted lines are shown for quark matter with 
a s = 0.5 and a s = 0.6 with bag pressure B = (155 MeV) 4 respectively . A remarkable feature of the state of affair 
is that there exist transition points for nuclear matter to quark at diffrent pressures and chemical potentials . These 
trasition points (P cr u, Aw) are at a s = 0.6 with (150 MeV/fm 3 , 1280 MeV) and at a s = 0.5 with (260 MeV/fm 3 , 
1445 MeV) showing dependance on a s and these also indicate the first order phase transition from nuclear matter to 
quark matter at differnt thermodynamical conditions . We also note that the phase transition seems to occur around 
the number density of about 5 times the nuclear matter density. These points also change under differnt bag pressure. 
Figure 3 shows the phase diagram with a different bag pressure B = (150 MeV) 4 . Here we found that transition 
point decreases with decrease of bag pressure, whereas the transition point shifts to a higher value with an increase 
in coupling constant a s .The early phase transition from nuclear matter to quark matter obviously implies that the 
interior of "neutron star" will usually consists of quark matter. We investigate this possibility in the next section. 
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FIG. 2. Pressure (P) versus chemical poten- FIG. 3. Same as Fig.2 with a different bag 
tial (fj.) for nuclear matter and for quark matter pressure (B). 
with various a s at constant bag pressure (B). 



IV. HYBRID STARS 



For the description of neutron star, which is highly concentrated matter so that the metric of space-time geometry 
is curved and one has to apply Einstein's general theory of relativity. The space-time geometry of a spherical neutron 
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star described by a metric in Schwarzschild coordinates has the form EQ] 

ds 2 = -e v{r) dt 2 + [1 - 2M(r)/r]- 1 dr 2 + r 2 [d0 2 + sin 2 Qd4> 2 ] (4.1) 
The equations which determine the star structure and the geometry are, in dimensionless forms |20|]2l| l 

dP(fr ) = g [e(fro) + P(rr )][M(rr ) + 47rqf 3 P(rr )] 



f/f 



f 2 [l - 2GM(fr )/r] 



M(fro) = 4™ / dr'f' e(r'ro), 



and the metric function, i/(r) is given by 

du(fr ) 



df 



2G 



[M(fr Q ) + Airaf 3 P(fr )} 
f 2 [l-2GM(fr )/f] 



In equations (4.2) the following substitutions have been made. 



e = e/e c , P = P/e c , f = r/r , M = M/Mq, 
where, with fi = 197.327 MeV fm and r = 3 x 10 19 fm, we have 



e c rg/M , G 



G/iM e 

To 



(4.2b) 



(4.2c) 



(4.3a) 



(4.3b) 



In the above, quantities with hats are dimensionless. G in equation (4.3b) denotes the gravitational constant with 
G = 6.707934 x 10" 45 MeV" 2 . 



In order to construct a stellar model, one has to integrate equations ( |4.2aD to (|4.2c|) from the star's center at r = 
with a given central energy density e c as input until the pressure P(r) at the surface vanishes. As stated in the last 
section, with any reasonable central density, we expect that at the center w e sh all have only quark matter . Hence we 
shall be using here the equation of state for quark matter through equation (3.4) with P(0) = P(e c )- We then integrate 
the TOV equations until the pressure and density decrease to their critical values at radius r = r c . For r > r c , we 
shall have equation of state for neutron matter where pressure will change continuously but the energy density will 
have a discontinuity at r — r c . TOV equations with equation of state for neutron matter shall be continued until the 
pressure vanishes. This will complete the calculations for stellar model for hybrid "neutron" star, whose mass and 
radius can be calculated for different central densities. 

In Fig. 4a, we plotted the mass of the star as a function of central energy density to examine the stability of such 
stars. As may be seen from the figure, dM/de c starts becomeing negative around 1500 MeV/fm 3 after which it 
becomes unstable and may collapse into black holes 0,^2| with the Chandrasekhar limit as 2.3M©. This yields 
stable hybrid star of mass M ~ 2.3-Mq with radius R ~ 13.5 km with a quark core around 4.2 km as seen from Fig. 4b. 
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Fig.4a 



Fig.4b 
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FIG. 4a. The mass of the hybrid neutron star FIG. 4b. The mass ratio (M/Mq) versus ra- 
(M/Mq) as a function of central energy den- dius (R) of the hybrid stars, 
sity (e c ). 



We also calculate the surface gravitational red shift Z s of photons which is given by |2C 

1 



Z s 



- 2GM/R] 



1. 



(4.4) 



In Fig 5 we plot Z s as a function of M/Mq. It is however possible that the discrete slowing down of pulsars due to 
the presence of two states of matter with various mass throwing some light on the above structure. Our graph shows 
there is a discontinuty of Z s around M/Mq — 0.4 indicating a peculiar behaviour of redshift with Mass of the hybrid 
star. 

Since the stars rotate about a centre, the relation between relativistic Kepler frequency and Newtonion Kepler 
frequency is given by 



Q k ~ 0.65ft c 



where 



Q c = x/M/R 3 



= 3.7 x 10 d 



' M/Mq 
(R/km) 3 



(4.5a) 

(4.5b) 
(4.5c) 



fifc is newtonian Kepler frequency balancing gravity with centrifugal force . The factor 0.65 is emparical and 
approximate. The figures Fig. 6 shows that 
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Fig.5 



Fig.6 
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FIG. 5. The surface gravitational redshift FIG. 6. The mass ratio (M/M Q ) versus ke- 
(Z 3 ) as a function of star mass (M/Mq) plerian frequency (Q c ) showing the mass ratio 

restricted to 2.3 

the frequency is higher for lower mass and decreases for higher mass i.e. M/Mq is about 2.3 . It clearly shows that 
we can not have mass of the hybrid star more than about 2.3 times that of sun . That also indicates (e.g. Fig. 4b) the 
radius of the star can not increase beyond 13.25 km . 



V. SUMMARY AND CONCLUSIONS 



We considered the equation of states taking into account the self-coupling interactions of a and u>- mesons. The 
inclusion of the quartic term to the uj meson field gives a soft equation of state. In our calculations, we used the 
TM1 parameter set, which has a capability to reproduce the known results of finite nuclei as well as of normal nuclear 
matter. Here also the TM1 parameter set gives a phase transition for hadronic matter and quark matter. The same 
parameter set predict the Chandrasekhar limit for Hybrid stars to be 2.3 Mq. In our calculation, we predict the inner 
quark core having a radius of about 4.2 kilometers, whereas the total radius of the hybrid neutron star is found to be 
13.5 kilometers as compared to the earlier result where Chandrasekhar limit is 1.58M Q and radius around 10 km JTl| , 
where the nuclear matter equation of state was calculated through the dressing of pion pairs. This is due to the 
contribution from a — ui mesons. One also notes that the redshift has discontinuity around M/Mq = 0.4, a peculiarity 
of hybrid stars. It is also observed that the Newtonian Kepler frequency of the hybrid stars can not increase beyond 
M/Mq = 2.3, showing a decrease with increase in M/Mq. Pulsars are expected to be stars of this type but the gross 
properties appear to be similar to what we believe regarding neutron stars. 
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